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We report on recent results obtained for the 3-loop heavy flavor Wilson coefficients in deep- 
inelastic scattering (DIS) at general values of the Mellin variable at larger scales of Q^. These 
concern contributions to the gluonic ladder-topologies, the transition matrix elements in the vari- 
able flavor scheme of 0{nfTp) and 0{Tp), and first results on higher 3-loop topologies. The 
knowledge of the heavy flavor Wilson coefficients at 3-loop order is of importance to extract the 
parton distribution functions and (M|) in complete NNLO QCD analyses of the world precision 
data on the structure function F2{x,Q^). 
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The massive Wilson coefficients in deep-inelastic scattering are known to be expressible in 
the limit of high virtualities » as convolutions between massive operator matrix elements 
(OMEs) and massless Wilson coefficients [1]. Here m denotes the heavy quark mass. The gen- 
eral structure of the Wilson coefficients to 0{a^) has been derived in [2]. These massive Wilson 
coefficients are in turn convoluted with parton distribution functions to obtain the heavy flavor con- 
tributions to DIS structure functions at leading twist. They have been calculated for the twist-2 
heavy flavor contributions to the unpolarized structure functions at leading [3] and next-to-leading 
order [4] ^ for general values of Q^. Since the massless Wilson coefficients are known by now at 
3-loop order [6], it remains to compute the OMEs analytically ?itO{a^), in order to obtain the mas- 
sive Wilson coefficients at NNLO. These coefficients will allow for a consistent NNLO analysis of 
the deep-inelastic world data at 2^>20GeV^, cf. [7]. 

In these proceedings, we discuss recent progress obtained in this direction. Our aim is to calcu- 
late all contributing OMEs for general values of the Mellin variable N. An important previous step 
towards this goal was the computation of the moments of the massive OMEs for = 2 . . . 10(14) 
contributing in the fixed and variable ^ flavor schemes [2]. The 3-loop heavy flavor corrections to 
Fi{x, Q^) in the asymptotic case were calculated in [9]. First results for general values of N for the 
color factor factors TpCA^F were calculated in [10] for two heavy quark lines of the same mass. The 
case of two different quark masses was considered in [10, 11] for fixed moments. Results for the 
color factors nfTpCA,F for general N were obtained in [12, 13] and the calculation of 3-loop ladder 
topologies was performed in [14]. Two-loop results up to 0{e) were obtained in [15]. Here the 
massive OMEs are computed for on-shell external massless partons. The case of a massive on-shell 
external fermion line was studied at two loops in [16] in case of QED. 

In the following we will describe the methods used to perform these computations. We gen- 
erate the Feynman diagrams using QGRAF [17]. After the numerators of these diagrams are con- 
tracted with appropriate projectors we end up with a large set of scalar integrals. Many of these 
integrals are calculated using a variety of approaches, namely, 

1. Modern summation algorithms, implemented in the Mathematica package Sigma [18]. 

2. The method of hyperlogarithms for convergent integrals, generalizing the method developed 
in [19] to one additional variable x. 

3. Mellin-Barnes integral representations [20]. 

4. The use of integration by parts identities [21] to express all integrals in terms of a small set 
of masters integrals. 

We will focus here on the first two methods and show a few examples. The Feynman diagrams 
with operator insertions may be turned into nested sums [22] . These infinite and finite sums may 
be solved using Sigma whenever they have a representation in terms of elements of difference- 
and product fields. This includes divergent diagrams, since the different poles and powers in e may 
be separated. Let us consider the scalar integrals associated with the ladder diagrams like the one 
shown in Fig. [1} In this diagram, the loop fermion is massive and the momentum of the external 

'For a precise implementation in Mellin space see [5]. 

^In using variable flavor schemes a correct scale matching is of importance [8]. 
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gluons is p, with = 0. We consider the case where all powers of propagators are equal to one, 
and in the numerator of the integral we only have the operator insertion (A • l)^ . The result after 
Feynman parameterization and calculation of the loop-momentum integrals turns out to be [14] 



i{A.pfa^,Sl f 



J _ -v 1' / -■s-'E j 

i^n — : — 3 l\a : 



(m2) 



2A2-4e 



where is the spherical factor Se = exp [f (/£ — ln(47r))] , and 

/ 3 , , 0(1-Wl-W2)w7''/V2''^^l-Wl-W2) 



1 +^1 



2-3e/2 



(1) 



1— W3 ^ 1— W4 
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Expanding the polynomial that appears raised to the Mh power, one can see that the w\- and 
W2-integrals, can be written in terms of an Appell hypergeomteric function. After an appropri- 
ate analytic continuation, we end up with the following representation of the integral in terms of 
multiple sums. 
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(3) 



We can now expand in £, and the resulting multiple sums can then be performed using the package 
Sigma. The result for this and other integrals can be written in terms of harmonic sums Sa [23] 
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and their generalizations S3{^;N) [24,25] ^: 



k=\ 



5,,2(t],^;A^) = 50 = 1, T],^e: 



k=i 



(4) 



Omitting the explicit dependence of the harmonic sums on N, we obtain 



4(A^+l)5i+4 
'(A^+l)2(A^ + 2) 
4(A^+l)5i -AN 

nT\ 

7A^ + 11 



C3 + 



252, 



1,1 



+ 



1 



(A^ + 2)(A^ + 3) (A^+l)(A^ + 2)(A^ + 3) 



c4 

-2(3A^ + 5)53,1-^ 



52,1+2 

5A^ 



(2A^ + 3)5i + 



I C _ ^ c2 

■(A^+l)(A^ + 2) A^+1 ^ 2 \ 
I, , 2A^ + 3 



5A^ + 6' 



53 + 



2{3N + 5)S'l 9 + AN 2 



{N+\){N + 2) 



+ ■ 



^2 H 5? + ^ — TTT, — ^ — zS\ 

N+\ ^ {N+lY{N + l) 



(5) 



This result was checked using MAT AD [27] for the fixed moments N = \ . . .10. Other, more in- 
volved, integrals calculated in a similar way were given in Ref. [14]. 

The second method we have used to compute the integrals is based on an algorithm originally 
proposed in [19]. It is applicable when the integral turns out to be finite, even in case for local 
operator insertions for a fixed integer value of the Mellin variable N. We have generalized this 
method to the case allowing for one non-vanishing fermion mass and local operator insertions 
in order to find the general A^-representations for convergent 3 -loop topologies. We work in the 
a-representation and obtain integrals of the form 

/4(A'^) ^ j "' j "^^2 da^ da4 das da^ da-j da^ ij2y2 ^ "IL'^'^ ■ 

The corresponding graph polynomials of a graph G are given by 

• U = Y,T ri/^r OJ/' where T denotes the spanning trees of G. 

• Dodgson polynomials [28] T arise from the operator insertions. The form of these poly- 
nomials will depend on the specific operator insertion we are considering. 

The integrals given by are projective integrals, where one a-parameter may be set to one 
eliminating the 5-distribution. The operators sit on on-shell diagrams which obey specific symme- 
tries. These are generally not obeyed by the operator insertion. The Feynman parameter integrals 
are now performed in terms of hyperlogarithms [19] L(>^,z) : C\r — )• C, where 

• r = {ao, CJi , aAf} are distinct points in C which may contain integration variables. 



^Cyclotomic and generalized cyclotomic harmonic sums and polylogarithms and tlieir relations have been treated 
in [26]. 



4 



New Results on the 3-Loop Heavy Flavor Wilson Coefficients . 



Abilio De Freitas 



• is a word over the alphabet 2t = {ao,ai ,...,0^}, where each letter a,- corresponds to a 
point a, . 

L( v^,z) is uniquely defined by the following properties : 



L({},z) = l, and L{0",z, 



1 



log"(z) for « > 1 



^L{{ai^},z) = -^L{^,z) for z £ C\L 
oz Z- Oi 

If is not of the form w = (0,0, • • • ,0), thenlimL(M^,z) = 0. 



(7) 



For example, L {ai,z) = log(z - a,) - log(a;). 

The hyperlogarithms satisfy shuffle relations, e.g. 

L{{ai,a2},z)L{{a3},z) = L{{ai,ai,a2},z) + L{{aua3,a2},z) + L{{ai,a2,ai,},z) ■ (8) 

The points to which the indices a, correspond may contain further integration variables. Using 
these properties after partial fractioning and integration by parts, one can express any primitive for 
expressions consisting of rational and hyperlogarithmic functions in terms of different hyperloga- 
rithmic functions. These primitives have to be evaluated at the respective integration limits. Due 
to the operator-insertions leading to power-type functions, the integrals do not fit directly into the 
framework of the algorithm for general values of N. In order to obtain the corresponding extension 
a generating function is constructed by the mapping. 



1 



k=0 



1 -xp{ai,--- ,an) ' 



(9) 




Figure 2: A 3-loop Benz diagram. 

Performing the Feynman-parameter integrations then leads to an expression which contains hyper- 
logarithms L„, in the variable x. Using this method, the scalar integral with all powers of propagators 
equal to one associated with the diagram shown in Fig. ^ corresponding to a Benz-type topology, 
yields 



/(x) 



1 



(^3 2L({-l},x) -2(-l+2x)L({l},x) -4L({l,l},x) 



;i+A^)(2 + A^)x 

-3L({-1, 0,0, l},x) +2L({-1, 0, 1, l},x) -2xL({0,0,l,l},x) +3xL({0, 1,0,1}, x) 
-xL({0,l,l,l},x) + (-3 + 2x)L({l, 0,0,1}, x)+2xL({l, 0,1,1}, x) -L({1,0, 1, 1, l},x) 



5 



New Results on the 3-Loop Heavy Flavor Wilson Coefficients . 



Abilio De Freitas 



-(5x-l)L({l,l,0,l},x)+xL({l,l,l,l},;c)-2L({l,0,0,l,l},x)+3L({l,0,l,0,l},;c) 
+2L({l,l,0,0,l},x)+2L({l,l,0,l,l},x)-5L({l,l,l,0,l},x)+L({l,l,l,l,l},x)| . 

(10) 

Finally, the A'^th coefficient of this expression in x has to be extracted analytically in order to undo 
the mapping This is achieved using the GetMoment s option of the package Harmonic Sums 
[25]. One may also use guessing-methods to obtain the corresponding difference equation based 
on a huge number of moments, cf. [29]. For a more complicated graph with non-trivial argument 
structure in x we were able to produce ~ 1500 moments [30]. One obtains from Eq. (10) 

1 J 648 + 1512A^+1458A^2_^744A^3_^212A^'^ + 32A^5_^2A^6 

- (A^+1)(A^ + 2)(A^ + 3)| (1+A^)3(2 + A^)3(3+A^)3 

2(-l + (-l)^ + A^+(-l)^A^) .3 , 1 e4 1, 

{l + 22N+mN^) 19 2 l+4N + 2N^ i ^ ^ (-9 + 4jV) 
2(1+A^)2(2 + A^) ^ 8 2 2(1+A^)2(2 + A^) ' 4 ^ 3(1 +A^) ^ 
^, ^.N^ {-l + 6N)^ 54 + 207N + 246N^ + \30N^ + 32N'^ + 3N\ 
5-2,1 + ^yTa^^^'1 + (i+Ar)3(2+A^)2(3+A,)2 

(-2 + 1N) 13 1 

+4^3^! - ^^^^ ^jv) '^2^1 + y 53^1 -7^2,iSi -753,1 + 1052,1,1 >■ (H) 

Another example, where this technique has been applied, is shown in Fig. ^. 




Figure 3: A second example of a 3-loop Benz topology. 



In this case the result is 



1 f2(l-13(-l)^ + (-l)^23+^+A^-7(-l)^A^ + 3(-l)^2i+^A^) 



1 . I (-If .3 (-l)^(3 + 2jV) ^ 5(-l)^ ^2 I 2(-l)^(3+jV) , 
{2 + Ny^ 2{2 + Ny^ 2{l+Nf{2 + N)^^ 2 ^ + (1 +A^)(2 + A^) ' 



I (-If (3 + 2jV) (-If, ,2 , 3(-l)^(4 + 3jV) , 2 ^ 

-12(-l)%C3 + ^y_^'^^^(^^^^^ ^ 
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4((-l)A'22+^-3(-2)^A^ + 3(-lf2i+^A^) /I , 

-•^1,2 T ) M — Ij '^2,1,1 



(1+A^)(2 + A^) ^'"V2 



, 2{-{-lf2^+^-l3{-2fN + 5{-lf2^+^N)^ fl 



+2(-»"f3^. P) + (i+A,)(2+A.) ir ' 

-2(-l)%,i,2(2.^,l)-(-l)%,l.l,i(2.i,l,lU. (12) 

Notice the presence of generalized harmonic sums and highly divergent factors in the limit 
— )■ oo, such as 2^. It can be shown, however, that the complete expression is convergent in this 
limit and possesses a well-defined asymptotic expansion for N ^ oo. In general neither the repre- 
sentation in individual nested sums or by iterated integrals shows this property, but a corresponding 
combination of terms does. 

We calculated the contributions of 0{nfTpCA^F) to all massive OMEs completely [12, 13]. 
Furthermore, first systematic results were obtained for the case of graphs containing two massive 
fermion lines with mi = m2. A typical graph is shown in Fig. H in the gluonic case. 




Figure 4: A 3-loop graph containing two massive fermion lines mi — OT2 and an operator insertion. 
One obtains 



1 



105£2 £ 



74A^3 _ 455^2^ 33 210 i 
44100(A^- 1)A^(A^+1) 210 



8903A^3 ^ 39537A^2 _ ^ i444oa^ + 36576 
2822400(A^+1)(2A^-3)(2A^- 1) ' 



^ 148176000(A^-1)2A^2(^+1)2(2A^_3)(2A^_1) +'^2(A^) + 3C2 

2-^N-\N-\)N{5N-6) (2N\ ( _f ^ , flMl^ 

3(2^ - 3)(2^ - 1) \n)\ ^, (2/) /3 + -e, ei) f j ' 



7=1 



(13) 



3 



Pi =1795487A^^ -7087789A^^ + 10654130//'^ -5797102//^ +6828839A^'^ - 16594069A^ 

+ 9651 144Af2 + 902160A^ - 1058400 . (14) 

Integrals of this type usually contain finite binomial and inverse binomial sums, which even may 
be nested. 

In conclusion, we have seen that the methods shown here allow us to obtain analytic expres- 
sions at general values of for 3-loop integrals contributing to the massive OMEs which could not 
be obtained by other methods before. We continue working on the set of integrals that we need in 
order to obtain aU necessary operator matrix elements. In particular, we are studying the possibility 
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also to extend the method of hyperlogarithms. Application of these methods to the more compli- 
cated case of non-planar integrals are underway. The package Sigma and related packages are 
continuously being upgraded to be able to meet the challenges that keep arising in this endeavor. 
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